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We propose a method to generate high field terahertz (THz) radiation with peak strength of GV/cm level 
in the THz frequency gap range 1-10 THz using a relativistic laser interaction with a gaseous plasma target. 
Due to the effect of local pump depletion, an initially Gaussian laser pulse undergoes leading edge erosion and 
eventually evolves to a state with leading edge being step function. Interacting with such a pulse, electrons 
gain transverse residual momentum and excite net transverse currents modulated by the relativistic plasma 
frequency. These currents give rise to the low frequency THz emission. We demonstrate this process with 
one and two dimensional particle-in-cell simulations. 


I. INTRODUCTION 

High field terahertz (THz) pulses with peak electric 
field strength > 1 MV/cm in the previously hardly ac¬ 
cessible THz frequency gap range 1-10 THz have gained 
wide scientific interest, since they offer attractive oppor¬ 
tunities for studies in condensed matter physics, biology, 
chemistry and photonics. As a unique and versatile tool, 
high field THz sources not only allow probe and observe 
the state of matter with high sensitivity, but also open a 
door to manipulate and control the electronic, ionic and 
spin degrees of freedom of matter both resonantly and 
nonresonantly. ^ 

Gurrently, high field THz pulses can be generated from 
large-scale accelerators using ultrashort electron bunches. 
For example, peak field of 44 MV/cm has been obtained 
via coherent transition radiation in the Linac Goher- 
ent Light Source (LGLS).^ In the laser-based schemes, 
THz pulses with field strength > 8 MV/cm at 1 KHz 
repetition rate have been generated via two-color laser 
filament at ion. ^ Single-cycle THz pulse with field strength 
up to 36 MV/cm can be generated by optical rectifi¬ 
cation of a midinfrared laser in a large-size nonlinear 
organic crystals assembly.^ In addition, other methods 
that have produced high field THz pulses include differ¬ 
ence frequency mixing process of two near-infrared lasers 
in second-order nonlinear crystals, but only with compo¬ 
nents > 20 THz,^ and optical rectification in lithium nio- 
bate (LiNbOs) crystals with tilted laser front, but mostly 
limited to < 1.5 THz.^ Besides, high field THz source can 
be generated from relativistic laser irradiated plasmas via 
various mechanisms. 

In our previous work, we have shown numerically that 
THz pulses with field strength > GV/cm can be gener¬ 
ated in plasma wakefieds by use of temporally tailored 
laser pulses.The mechanism is related to the excita¬ 
tion of transverse plasma current via electron residual 
momentum left by the temporally tailored laser pulse. 
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For an electron initially at rest, the transverse momen¬ 
tum pz normalized by rrieC after the interaction with a 
2 :-polarized laser pulse can be obtained as an integral of 
the transverse field Ez'^‘^ 
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where az is the normalized vector potential, r] ^t — xjc^ e 
is the elementary charge, rrie is the electron mass, and c is 
the speed of light in vacuum. After the interaction with 
a normal Gaussian pulse, the above integral is zero, and 
hence electrons gain no net energy. However, if the laser 
pulse has such a steepened edge that it becomes a step 
function, a large transverse residual momentum as 
well as vector potential will be generated. Gonsidering 
the one-dimensional (ID) wave equation^^ 


= (^) Szix,t) (2) 

with the source term 
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where aj is the vector potential of the THz wave, 
^po = (47rnoe^/me)^/^ is the background plasma fre¬ 
quency, ne{x,t) is the plasma density normalized by the 
initial density no, and 7 is the relativistic Lorentz factor. 
The normalized THz electric field by can then 

be obtained as 
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where Ig = cjujpQ and ujq is the laser angular frequency. 

In this work, we demonstrate that THz emission from 
laser plasma wakefields can be self-consistently realized 
with a normal Gaussian laser pulse, unlike previously 
the laser pulse with a step-function edge is introduced 
artificially.^^ The way to achieve the required pulse shape 
is to utilize the effect of localized etching. When a high 
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FIG. 1. (Color online) Spatial profiles of the electron density 
Ue (green) and transverse field Ez (red) at (a) ti = 80To and 
(c) t 2 = 3 I 8 T 0 . (b) and (d) are the corresponding spatial 

distribution of the transverse momentum Pz for (a) and (c), 
respectively. The inset of (c) shows a close-up for the laser 
pulse shape. 


intensity laser pulse propagates in under dense plasmas, 
the leading edge of the laser pulse can get eroded and 
depleted via the local pump depletion mechanism, and 
thus eventually evolves to a state with a step-function 
front. The net transverse current exited by such a pulse 
in the plasma wakefields, modulated by the relativistic 
plasma frequency, should be capable of generating low 
frequency THz pulses. 


II. SIMULATION RESULTS AND DISCUSSIONS 

To illustrate the THz generation process, ID particle- 
in-cell (PIC) simulations have been carried out using the 
Virtual Laser Plasma Laboratory (VLPL) code.^^ A cir¬ 
cularly polarized laser with a Gaussian temporal profile 
a{t) = ao exp(—) is normally incident along the x- 
axis, where ao and tq are respectively the laser normal¬ 
ized peak amplitude and the pulse duration. We note 
that linearly polarized laser also works. Here, we take 
ao = eEoImeCUJo = 30 with Eq the laser electric field am¬ 
plitude. The pulse duration is about 11 fs full-width at 
half-maximum (FWHM). The laser wavelength is taken 
to be Ao = 0 . 8 /im. A fully ionized plasma is initially 
located between 50Ao and 350Ao- For the laser intensi¬ 
ties and plasma densities we considered here, the usual 
gas jets routinely used in experiments can be employed, 
such as hydrogen, helium and nitrogen gases.For 
simplicity, a homogeneous plasma slab with a density of 
no = O.OlUc is used, where ric = ^ 10 ^^ 

W/cm^ is the critical plasma density. We also checked 
inhomogeneous density profiles with density gradients 
and find this mechanism also works. The simulation box 
length is 550 Ao- The grid size is Ao/lOO with each cell 
filled with 10 macroparticles. 

Figure 1 shows the snapshots of distribution of the 
laser pulse E^ and the electron density as well as 
the transverse momentum of electron Pz at two different 
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observed in front of the right plasma-vacuum bound- 


times: ti = 8 OT 0 and t 2 = 3 I 8 T 0 , where Tg is the laser 
period. In Fig. 1(a), an electron density peak is formed 
at the front of the laser pulse by snow plowing the elec¬ 
trons forward due to the laser ponderomotive force. The 
peak electron density is larger than the initial background 
density, but still much smaller than the critical density. 
The laser pulse almost holds its initial symmetric tem¬ 
poral shape. The corresponding transverse momentum 
Pz in Fig. 1(b) only exits inside the laser pulse. Behind 
the laser pulse, one finds Pz ~ 0. From Fig. 1(c), how¬ 
ever, one can see that at a later time, the density peak 
at the front of the pulse becomes narrow and sharp. The 
amplitude of the density spike is above the critical den¬ 
sity. This results in an efficient erosion of the leading 
edge of the pulse due to localized energy depletion mech¬ 
anism within the narrow density spike region. As the 
laser pulse propagates, the density spike is continuously 
pushed by the laser ponderomotice force in front of the 
pulse. Consequently, the laser front coinsides with the 
density spike where it depletes its pump energy. Due to 
the etching of the front, the laser pulse evolves to a state 
with an ultra-steep front which is approximately a step 
function (see a close-up for the laser pulse shape in the 
inset of Fig. 1 (c)). As the laser front continuously etches 
backwards, the sharp shape of the leading edge is kept. 
Although there are some other effects like group velocity 
steepening that can lead to the formation of steep pulse 
front,these processes are not so efficient. With this 
pulse shape of step-function front, the integral in Eq. (1) 
is no longer zero. As a result, the electrons acquire a net 
transverse momentum after the interaction with the laser 
field. This residual momentum left behind the laser pulse 
can be clearly seen in Fig. 1(d). Due to the residual mo¬ 
mentum, the plasma electrons keep on a free oscillation 
with the background plasma frequency Then a low 
frequency transverse net current is built up in the plasma 
wakefields and emits both forward and backward low 
frequency electromagnetic radiation through the plasma 
slab.^^ We noth that the second density peak at around 
260 Ag behind the laser pulse in Fig. 1(c) is a trapped 
electron sheet injected into the wakefield. This electron 
sheet may generate an ultrafast XUV pulse in the pres¬ 
ence of the residual transverse momentum. 
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FIG. 3. (Color online) (a) x — pz and (b) x — px phase 
space distributions of the electrons, (c) Spatial distribution of 
the relativistic 7 -factor, (d) Spatial profiles of the relativistic 
plasma frequency ojprei (blue) and the transverse electric field 
Ez (red) in the skin layer. The black dashed line marks the 
initial plasma frequency The time is t = 360To. 


FIG. 4. (Color online) (a)-(b): Peak THz field strength as 
a function of (a) the driven laser amplitude ao and (b) the 
initial plasma density no. (c)-(d): THz spectrum dependence 
on (c) ao and (d) no. The plasma density is no = O.Olnc for 
panels (a) and (c). The laser amplitude is ao = 30 for panels 
(b) and (d). 


The temporal waveform and the corresponding Fourier 
spectrum of the low frequency emission observed in the 
forward direction are shown in Figs. 2(a) and (b), respec¬ 
tively. A single-cycle THz pulse is obtained, which lasts 
for a few hundred femtoseconds. The peak field strength 
of the THz pulse reaches about 24 GV/cm. The spectral 
shape shows that most of the pulse energy concentrates 
in the desired THz frequency gap range of 1-10 THz, with 
a peak frequency of about 3 THz. 

The initial plasma frequency is ~37.5 THz, which 
is greater than the central frequency of the observed THz 
emission. This can be attributed to the change of the 
plasma frequency due to relativistic effects. As a re¬ 
sult of the interaction with the relativistic laser, elec¬ 
trons reach nea r light spe ed and gain a large Lorentz 
factor (7 = 1 / — [v/cY). Therefore, the mass of elec¬ 
tron is increased and the effective plasma frequency is 
reduced. The relativistically corrected expression for the 
plasma frequency is cjprei = (47re^ne(x)/7me)^/^. We 
plot the electrons’ 7 -factor along with the phase space 
distribution at time t = 360To in Figs. 3 (a)-(c). The 
7 -factor is obtained by averaging the values in each cell, 
i.e., 7 =< 7 >ceii- In the skin layer near the plasma 
rear surface, large residual transverse momentum can be 
seen, which can excite net transverse currents radiating 
electromagnetic pulses (see Fig. 3(a)). In this plasma re¬ 
gion, the Lorentz factor 7 of the bulk electrons reaches up 
to 100 (see Fig. 3(c)). In addition to the bulk electrons, 
there are some high-energy electrons with 7 > 100 , which 
can also be seen in Fig. 3(b). These electrons are re¬ 
fluxing backwards in the electrostatic field due to charge 
separation with the ions. Unlike the bulk electrons, they 


contribute little to the plasma frequency. Figure 3(d) 
shows the profiles of the relativistic plasma frequency 
<^prei(^U) fhe transverse electric field Ez{x^t) in the 
skin layer. The initial undisturbed plasma frequency ujpo 
is marked as the black dashed line. Due to the relativistic 
effect, the effective plasma frequency is largely decreased 
to be around a few THz, which is in the THz frequency 
gap range. The net transverse current modulated by this 
plasma frequency leads to the emission at the low THz 
frequencies. 

We find the present scheme works for a wide range 
of laser and plasma parameters. Figure 4 presents the 
results of a series of simulations, where we take no = 
O.Olnc in panels (a) and (c), and ao = 30 in panels (b) 
and (d). Figs. 4 (a) and (b) shows the peak amplitude 
of the THz emission increases with increase of the driver 
strength ao and the plasma density no- Higher ao and no 
both lead to a higher density spike in front of the laser 
pulse, which results in a more efficient local depletion 
of the laser energy, and thus a stronger radiation field. 
A rough scaling of the THz strength dependence on the 
laser plasma parameters can be obtained using the above 
ID model. Assuming the source is sharply localized and 
exponentially decays in the skin depth /g, the integral 
of Eq. (4) can be approximated by ej ujpSz/2ujQ}^'^^ 
Concerning the source term, we further assume that the 
driver laser is nonevolving. With the assumption of the 
quasistatic approximation,^^ the plasma fluid quantities 
can then be expressed as ng/no = [ 7 ^ + (l+(/)^)]/ 2 (l + (/))^ 
and 7 = [ 7 ^ + (1 + 0 ^)]/ 2 (l + 0), where 75 _ = 1 + a^ 
and 0 is the electrostatic wake potential.In the near 
wavebreaking nonlinear limit (1 + (/>) ^ 1 / 7 ^, where the 
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plasma Lorentz factor 7 ^ = (1 — v‘^1 with Vp the 

phase velocity of the plasma wave. Then the THz field 
can be obtained as ej oc ujpoa^jp. 

The THz spectrum dependence on the laser intensity 
and plasma density are shown in Figs. 4(c) and (d), re¬ 
spectively. With increasing the laser intensity, both the 
spectrum width and the central frequency decrease, and 
the peak spectral strength increases. The red shift of the 
central frequency is due to a lower relativistic plasma fre¬ 
quency, as a consequence of a larger 7 -factor. Since the 
emitted pulse is single cycle, a lower-frequency (longer- 
wavelength) pulse has a longer temporal duration, and 
thus a narrower bandwidth. From Fig. 4(d) one can see 
that the central THz frequency increases with increasing 
the plasma density, which is due to an increased plasma 
frequency. These results may offer a way to obtain a 
tunable THz source. 

In the following we present 2 D PIC simulation results 
to check the present scheme working in multi-dimensional 
geometries. The size of the simulation box is 550Ao x 6 OA 0 
in the xy plane and the grid step is Ao/50 x Ao/20. The 
laser amplitude is ao = 30 and the spot size is 15 jam 
FWHM. The plasma density is no = 0.02nc. Figures 
5(a) and (c) show the 2D snapshots of distribution of the 
transverse electric field and the electron density at 
time t = 3 OOT 0 , respectively. The corresponding line-out 
plots along T = 0 are shown in Fig. 5(b) and (b). As 
in the ID case, the laser pulse drives an electron density 
spike at the front. The leading edge is eroding and be¬ 
coming ultra-steep. Figures 5(e) and (f) show a snapshot 
of the electric field Ez, i.e., the distribution of the forward 
emission, and its corresponding Fourier spectrum along 
T = 0 at time t = 500To. The inset of Fig. 5(f) shows the 
line-out plot of the THz electric field along T = 0. One 
can see that the emission fields also have strong peaks in 
the THz frequency gap range of < 10 THz. The peak field 
strength is about 0.3 GV/cm. The reduction of the field 
strength can be attributed to multi-dimensional effects 
such as transverse electron spreading and wave diffrac¬ 
tion. In the ID cases, the dense electron spike in front 
of the laser pulse maintains for a long time, provided the 
laser ponderomotive force is much larger than the elec¬ 
trostatic restoring force. In the multi-dimensional cases, 
the transverse spreading effects reduce both the average 
density and lifetime of the electron spike, which conse¬ 
quently impair the efficiency of the laser front etching 
and THz generation. Further increasing the plasma den¬ 
sity and focusing the THz emission may result in a higher 
field strength. 
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FIG. 5. (Color online) Results from 2D simulations, (a) A 
snapshot of the transverse electric field eEz/me(Xoc distribu¬ 
tion and (b) the corresponding line-out plot along T = 0 at 
time t = 3 OOT 0 . (c) A snapshot of the density Ue distribution 
and (d) its line-out profile along Y = 0 at t = 300To. (e) A 
snapshot of the forward emission field eEzIrrieU^oc distribu¬ 
tion and (f) the Fourier spectrum of the field along T = 0 at 
t = 5 OOT 0 . The inset of panel (f) shows the line-out plot of 
the THz field along T = 0, which was smoothed to get rid of 
the numerical noise. 


local pump depletion effect, the laser front is gradually 
etched to become a step-function state. Plasma electrons 
gain a net residual momentum after interaction with such 
a pulse and excite transverse oscillations modulated by 
the relativistic plasma frequency, which emits the low 
frequency THz pulses. This result opens the way for ex¬ 
perimental approaches. This tunable THz source with a 
peak field strength reaching GV/cm level and spectrum 
mainly within the THz gap range, should be of interest 
for high-field THz science and applications. 
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